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We study the 3D topology of Rayleigh- Taylor (RT) and Richtmyer-Meshkov (RM) single-modes, 
which includes bubbles, jets and saddle points. We present an analytic description of the interface 
as a whole, for arbitrary time-dependant acceleration g{t). The dependance of morphology on the 
lattice - Hexagonal, square or triangular -of bubbles are investigated. RM accelerations in the case 
of a large density ratio produce jets well separated from each other while, in RT case, jets are 
connected by liquid sheets. We compare our analytic results to numerical simulations. 
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The RT and RM instabilities (RTI and RMI) play an 
important role in astrophysics in ion Q and laser 

1^ inertial confinement fusion (ICF), in shock-tube mix- 
ing 1^,^ and in chemical, nuclear or thermonuclear com- 
bustion [||. The asymmetry caused by RTI and RMI 
in spherical implosions strongly alters the neutron yield 
and energy gain in ICF targets [|,|. Recently, it was 
proposed to use experiments on very powerful existing 
and future laser systems such as Omega and the National 
Ignition Facility in the USA or the Laser MegaJoule in 
France for modeling the unstable explosion of supernovae 
(SN) or unstable expansion of SN-remnants [|l],|. 

Flow with two shocks and mixing between them (similar 
to SN-flow) takes place during an expansion of detona- 
tion products after an explosion [HJ|] (compare with [0 ) . 
RTI/RMI are also significant for other astrophysical ap- 
plications such as planetary nebulaes, Wolf-Rayet stars 
and magnetospheres of neutron stars 1^. The physical 
origins of exchange instabilities (RTI and RMI) are con- 
nected with baroclinic generation of vorticity j|,|ll[] . RTI 
is driven by buoyancy (see reviews |^,|l^). RMI occurs 
after the passage of a shock wave through surface corru- 
gations 1,1,00. 

Configurations of 3D single-mode perturbations can be 
represented by bubble lattices having various geometrical 
symmetries: hexagonal (B6), square (B4) or triangular 
(B3). The B stands for "bubble" and the digits 6, 4 and 
3 correspond to the number of bubbles adjacent to the 
chosen one. Saddles and jets also form lattices. In the 
lattice B6 a jet J has three neighbouring jets. Therefore 
the lattice B6 is the same time as the lattice J3. 

Our goal is to describe 3D phenomena and their de- 
pendence on both lattice symmetry and acceleration pro- 
file, g{t). This approach is needed to understand phe- 
nomena occuring in ICF, astrophysics, etc., because real 
unstable flows are three-dimensional and the accelera- 
tion, g(t), satisfies neither the RTI nor the RMI condi- 
tions |13[| . To begin our study, let us compare the ex- 
amples in Fig.|l| (2D) and Fig.| (3D). They represent, 
at timc|^ t — 1 , the position of the periodically per- 
turbed interface rj in the RM case. At t = 0, the 
interface ri{x,0) = t^oo (resp. rj{x,y,0) = rjoci) corre- 
sponds to the plane z = and the potential of the ve- 
locity perturbation is ip{x,z,t = 0) cx cosxe"^ [resp. 
(p{x,y,z,0) cx (cosx -f cos?/)e"'^] in 2D (resp. 3D/B4). 
These initial conditions are called "standard" because 
they are widely used for small initial perturbations. The 
heavy fiuid, with density p, is above (z > tj) the light fiuid 
with density pi . Jets of dense fluid move down while light 
bubbles rise up. Analytically we consider the case with 
density ratio p, = pi/ p —t Q. Our method of calculation is 



^ Space is normalized by wavenumber k and thy k and initial 
amplitude of velocity wo (RM) or k and g (RT) . 



presented below. 

In the 2D case (Fig. P [resp. 3D case (Fig. 1)], 
the points B and J (resp. B, S and J) are very impor- 
tant since they correspond to stagnation points. At these 
points, the velocity of the fluid becomes zero relative to 
the surface rj. 

FIG. 1. Rippled structure of 2D solutions. Periodic se- 
quence of parallel valleys. There is a chain ..-B-J-B-,]-.. of 
tips of bubbles B and jets J. 



FIG. 2. Geometry of a single 3D mode in a square lat- 
tice (B4). Periodic cell structure in the horizontal plane or 
two-dimensional crystal formed from points of tips of bubbles 
B, jets J and saddles S marked by circles. 

The 3D/B4 flow is invariant relative to QO^-rotations 
around B and J vertical axis and to 180°-rotations around 
S axis (Fig. I). 

For small density ratios, /i <C 1, the whole flow is 
clearly divided into two qualitatively very different parts. 
The first is the bubble envelope imprinting into the dense 
fluid. Bubbles brake the initially continuous dense fluid 
and produce jets. Above the envelope, the dense fluid is 
still in a contiguous state. The second part is the ejecta 
(jets) pinched and driven down by the imprinting bub- 
bles. The points B belong to the contiguous fluid and 
the points S and J to the ejecta. The 3D pattern of the 
ejecta is rather complicated (Fig. ||). The ejecta consists 
of 1) wall type jets going down from B to S and 2) leg 
type jets going from S to J. The bubble has the shape of 
a well in the dense fluid. This well transforms into walls 
or skirts around the points B. 

The relevant topological measure of the form of the 
ejecta is the geometrical ratio 

r(t) = Azbs/^zbj, with AzMN = ^m(0 - ZNit), (1) 

where zb,zs and zj are the vertical positions of the 
points (Fig. F is the ratio of the length of "skirt" 
to the length of "legs" . 

FIG. 3. Numerical results for RM (the triad in the upper 
panel, t — 20) and RT (lower panel, t — 9) cases for the 
lattices B6, B4 and B3 from left to right. The ratio is /i = 0.1. 

The six direct numerical simulations (DNS) in Fig. ^ 
were done by a grid-characteristics method [|l^. Several 
interesting works are devoted to DNS of 3D RTI/RMI 
flows 

Before presenting our results, let us give the bound- 
ary conditions and spectral decomposition. The motion 
is described by a velocity harmonic potential (p (the vor- 
ticity is concentrated at the interface). Classical [plp^] 
kinematic and dynamic conditions are 

r]t = w\ - rixu\ - riyv\, v ^ {u,v,w}, r/t = dT]/dt, (2) 
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ipt\=v'^\/2 + g{t)r], f\ = f\r, = f\z=v{x,v,t)- (3) 

We represent the potential by a Fourier series and the 
surface 77 near stagnation points by a Taylor series 



N 



(p{x,Z,t) = ^ipn{t)c 



-nAz 



(4) 



JV 

E 

n=l 



where Cnx = cosnx and N is the truncation number. It 
defines the order of approximation of conditions (Hj^). 

2D solution. The expansion (Q) satisfies Aip — 0. The 
expressions (^,^ introduce geometrical rjQ, Ki, .., Kn and 
velocity tfi, ■■, ^pn unknowns for the ordinary differential 
equation system in our method of asymptotic collocations 
(MAC). We say asymptotic collocations because of the 
close connection to the method of ordinary collocations 
in which boundary conditions are approximated in a set 
of discrete points {xi}. The equations of the MAC appear 
asymptotically when all points Xi tend to point B, or S, 
or J [plp^. In 2D, these equations, for < 6, were first 
derived and integrated in in case of bubbles. 

The Fig. |] presents the A^ = 5 solution. Here, for 
the first time, we use high-order MAC to study jets. 
We carefully describe acceleration of the jet and for the 
RM case obtain very accurate values for its terminal 
(asymptotic) velocity, w,j{t = oo,N ^ 5) = —1.923 
{wj{t) = ?7o = "nixj^t)), standard initial conditions. The 
accuracy of this value is e(l) = lO"^", e(2) = 10~^ '^, 
e(3) = 10-316 and e(4) = IQ-^is^ ^here e{N) = 
\wj{po,N + 1) — wj(oo, A^)|/|wj(oo, A^)|. Since the er- 
ror is very small, the MAC may be used to control the 
accuracy of other methods. 

3D solution. For the B4 lattice we have 



,(S4) 



(x, y, 2, <) = E E 



,(t)a;2"y2™/(2n)!(2m)!. 



where < n -I- to < A^, 



c-nm 



expl^-qimAz), Koo{t) = 
im? . The unknowns in the system 
are ?7o, A'nm, </5rim- For A^ = 1, the indices (n,m) are 10 
and 01. For A = 2, they are 10, 01, 20, 11 and 02. Points B 
and J (but not S) are symmetric, ipnm = <^mn and Knm = 
Kmn- At the lowest order A^ = 1 (Layzer approximation 
1^,^-21 1) the unknowns, at points B and J, are 770, A", w, 
where K = A'lo = ATqi is the curvature and w = — tpio — 
—(^01 is the velocity of a bubble or a jet. The system 
A = 1 for the B4 lattice, valid in points B and J, has 



been derived and solved for the RM case in |16|. The 2D 
A' = 1 system was considered in [p[p"^-|2l| and 3D A = 1 
systems for B6 and B4 cases were examined in |P,|l6|,p0i . 



This is the first time we use MAC to describe the 
dynamics of all three kinds of points (B, S and J) at 
order A = 2. It appears much more complicated than 
A^ = 1 (compare also with 2D, A^ < 5). The termi- 
nal velocities of RM jet and saddle are wj{oo,N = 
1) = -V2, wj{oo,N = 2) = -1.698 and ws{oo,N = 
1) = -0.572, ws(oo,iV = 2) = -0.512. This gives an 
accurate asymptotic value of the geometrical ratio (|^), 
r ws{oo)/wj{oo) as t ^ 00. Fig. |^ gives an example 
of second order solution. 

For the B6 and B3 lattices, we have in first order 



^(B6,3) ^ [<^+0(t)c+ + (^-0(<)c- + ip+-it)c]e- 



Az 



where c± = cosA:±f, c = cosPf, = {1/2, ±V3/2, 0}, 
fc'^ = {1, 0, 0} and = ip~'^ = (^^ for B and J points. 

The first order dynamical systems for points B and J 
for all three lattices B6, 4 and 3 are the same ||l6| : rja — w, 
w = + Ag{t)K]/2{l + 2K), K = -{1 + 4.K)w/2, 

where / = df/dt. In this system, we emphasize that g(t) 
is an arbitrary function. Eliminating t between first and 
last equations, we obtain dK/dr]Q = —1/2 — 2 AT. The 
solution is 



X(77o) = -l/4 + cxp(-2?7o)/4. 



(6) 



The 2D analog is dK/df]Q = -1 - 3K and K ^ -1/3 + 
exp(— 377o)/3. The solution (|^) has linear asymptotes for 
1 770 1 <C 1 and tends to —1/4 for I770I ^ 1. In the linear 
stage we have 77 « r]o{cix + ciy)/2 (B4), 77 « 770(0+ +0^ + 
c)/3 (B6,3) and K -770/2. 

In RM case, the differential system has integrals m& 



1 + 2K 



1 , V2 + 4K + VT+4K 
—= m ■ 



1 + 4.K y/2 (V2-h l)exp(-\/2) 



Wot, (7) 



^-1 + ^lnp^^^^U.ot. (8) 
w 2V2 \V2-l V2wo + wj 

Substituting solution AT (770) (|^) into (0) we obtain 



1 + e^'fo 1 Vl + e-2r,o p-no 

— '-71^" VfTi = "°'- 



In (|7|]9|) Wo is the initial velocity of B or J points. We 
will assume that this initial velocity for the tip of bubble 
B is equal to 1 for all lattices. Then initial velocities 
of points J are -1/2,-1 and —2 for lattices B6,4 and 3 
respectively. Similarly to ([7[]9|), the 2D solution has been 
obtained previously P, p^|jl9|]2l| ] . The integrals (0^) give 
rio{t), K{t), w{t) in analytic form for all times from initial 
to asymptotic state. They are valid in the RM case for 
B and J points of B6,4 and 3 lattices. 

It is very suprising that the relation AT (770) (|^) between 
the main geometrical characteristics 770 and K for a bub- 
ble penetrating into the dense fluid does not depend upon 
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g{t) for N = 1. This means that the relation is only 
weakly dependent on g(t) in the general case with arbi- 
trary N. 

The system for = 2 is rather long and can not be 
given here. However, for saddles and A'^ = 1, we have: 



Vo 



f3, 



K = -{l + 3K)a + K(3, 



-Qa + (1 + 3Q)/3, 



(1 + K)a + Kl3 = a^ + g{t)K, 



-Qa+il + Q)f3 = 13^ +g{t)Q 



(10) 



(11) 



(12) 



(13) 



(14) 



in the B4 case with ip — {—acix + (3ciy)e . For B6,3 
we have 



?7o = 2a - 7, 



2k = -(1 + QK)a + 2(1 + 2,K)-f, 



2Q = -{3+lOQ)a + 2Qj, 



(15) 



(16) 



(17) 



2(1 + AK)a + 4(1 + K)j = (a - 2jf + 4g{t)K, (18) 



2(3 + 4Q)d + 4Q7 = 9a^ + Ag{t)Q, 



(19) 



where a and 7 are the amplitudes of potential ip = 
[—a{c^ + c~) + 7c]e~^^ written with respect to the 
point S. In the standard case, initial conditions are 
a(0) = 7(0) - -1/3 (B6), a(0) = /3(0) = -1/2 (B4), 
a(0) = 7(0) = 2/3(B3) and ?7o(0) = K{0) = Q(0) = 
(B6,4,3). From ^ and (|l|), initial velocities of the 
saddles are -1/3 (B6), (B4) and 2/3 (B3) - to be com- 
pared with initial velocities of jets. Although, the sys- 
tems are the same for B6 and B3, the initial conditions 
differ. From systems (|lO-14) or (p^p^ we found numer- 
ically the trajectories of saddles ?7o(i) and the evolution 
of the curvatures K{t),Q{t). The terminal velocities of 
saddles in RM case are ws{oo,N = 1) = -0.748 (B6), 
wsioo, 1) = -0.572 (B4) and ws{oo, 1) = -0.196 (B3). 

In Fig. ^ we present the time variation of the geo- 
metrical ratio r. Initial values of the ratio are TBsit ~ 
0) = 8/9, rB4(0) = 1/2 and ri33(0) = 1/9. Simulations 
(/i — 0.1) fit rather well with the theory (/i 0) although 
/i differs. We observe that the agreement between theory 
and simulation is better for the RMI than for the RTI. 
The increase of N significantly improves the accuracy 
(for the B4 system, the N = 2 curve and the one coming 
from simulation are very close for the RMI) . The function 



g{t) influences therefore the evolution of this ratio. The 
morphology of the ejecta mainly depends on the type of 
lattice. The shortest skirts are obtained for the B3 lat- 
tice (lower set of curves) which produces powerful and 
fast jets (Fig. 3). Previously, it has been shown |l^ that 
for random (turbulent) cases, the patterns are similar to 
the B6 and B4 lattices (B3 lattices occur only for special 
conditions). Moreover, B3 type structures may appear 
after reshock and the corresponding rephasing because 
B6 bubbles transform into B3 jets. 

In summary, we have considered the effects of lattice 
and time-dependant acceleration on the evolution of the 
interface. We found that the shape of imprinting bub- 
bles is very weakly dependent on these factors, but the 
position of the bubbles depends on the history of the ac- 
celeration. At the same time, the shape and dynamics of 
the ejecta appear very sensitive to both factors. 
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RBRF (00-01-00250, 99-02-16666) and from Landau In- 
stitute - CNRS jumelage (N.A.I.) . 

FIG. 4. Effect of lattice symmetry and acceleration history 
on the evolution of the shape of rj. The upper, middle and 
lower sets of curves correspond respectively to the B6, B4 
and B3 lattices. The grey and black curves correspond re- 
spectively to numerical simulations - see Fig. 3 - and to the 
analytical approach with TV = 1. The dashed curves have 
been obtained from theory with N = 2. The curves labelled 
by 1 (resp. 2) describe the RTI (resp. RMI). 
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